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SUMMARY 

Most  published  work  relating  to  dynamically  loaded  journal  bearings  are 
directed  to  determining  the  minimum  film  thickness  from  the  predicted  journal 
trajectories.  These  do  not  give  any  information  about  the  subsynchronous  whirl 
stability  of  journal  bearing  systems  since  they  do  not  consider  the  equations 
of  motion.  It  is,  however,  necessary  to  know  whether  the  bearing  system  opera¬ 
tion  is  stable  or  not  under  such  an  operating  condition. 

The  purpose  of  the  present  paper  is  to  analyze  the  stability  characteris¬ 
tics  of  the  system.  A  linearized  perturbation  theory  about  the  equilibrium 
point  can  predict  the  threshold  of  stability;  however  it  does  not  indicate 
postwhirl  orbit  detail.  The  linearized  method  may  indicate  that  a  bearing  Is 
unstable  for  a  given  operating  condition  whereas  the  nonlinear  analysis  may 
indicate  that  It  forms  a  stable  limit  cycle.  For  this  reason,  a  nonlinear 
transient  analysis  of  a  rigid  rotor  supported  on  oil  journal  bearings  under 
(1)  a  unidirectional  constant  load,  (2)  a  unidirectional  periodic  load,  and 
(3)  variable  rotating  load  are  performed. 

In  this  paper,  the  hydrodynamic  forces  are  calculated  after  solving  the 
time-dependent  Reynolds  equation  by  a  finite  difference  method  with  a  succes¬ 
sive  overrelaxation  scheme.  Using  these  forces,  equations  of  motion  are  solved 
by  the  fourth-order  Runge-Kutta  method  to  predict  the  transient  behavior  of  the 
rotor.  With  the  aid  of  a  high-speed  digital  computer  and  graphics,  the  journal 
trajectories  are  obtained  for  several  different  operating  conditions. 


INTRODUCTION 

There  are  two  principal  approaches  by  means  of  which  one  can  analyze  the 
whirl  Instability  of  a  rotor  supported  on  fluid  film  bearings.  These  are: 

(1)  linearized  (perturbation)  method  and  (2)  nonlinear  transient  analysis.  In 
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the  linearized  method  a  small  perturbation  of  the  journal  center  about  Its 
equilibrium  point  is  given.  The  stiffness  and  damping  coefficients  are  deter¬ 
mined  after  solving  the  basic  differential  equation.  These  coefficients  are 
then  used  in  the  equations  of  motion  to  find  the  critical  mass  parameter  and 
whirl  ratio.  The  mass  parameter,  a  function  of  rotor  speed,  is  a  measure  of 
stability.  The  nonlinear  transient  analysis,  on  the  other  hand,  gives  the 
journal  locus  and  from  this  one  can  know  about  the  system  stability. 

Stability  analysis  of  finite  journal  bearings  by  the  linearized  method  has 
been  given  by  Allaire  (ref.  1),  whereas  Akers,  Mlchaelson,  and  Cameron  (ref.  2) 
studied  the  same  bearing  configurations  using  a  nonlinear  transient  approach. 

In  reference  2  It  was  shown  that  under  certain  operating  conditions  the  journal 
motion  was  bounded  and  could  form  a  limit  cycle. 

The  aim  of  the  present  paper  Is  to  study  theoretically  the  stability  char¬ 
acteristics  of  finite  oil  journal  bearings  under  dynamic  load  using  a  nonlinear 
transient  method.  A  few  papers  (refs.  3  to  5)  deal  with  the  dynamically  loaded 
bearings  to  predict  the  journal  locus.  As  these  do  not  consider  the  equations 
of  motion  for  the  prediction  of  the  position  of  journal  center,  they  cannot 
indicate  whether  the  bearing  system  Is  stable  or  not.  However,  these  are  use¬ 
ful  for  estimating  minimum  film  thickness  of  dynamically  loaded  bearings.  The 
dynamical  equations  of  motion  are  solved  by  fourth-order  Runge-Kutta  method  to 
find  eccentricity  ratio  and  attitude  angle  and  their  derivatives  for  the  next 
time  step.  These  values  are  then  Introduced  in  the  two-dimensional  time- 
dependent  Reynolds  equation  to  find  the  hydrodynamic  forces.  Following  the 
above  approach  a  nonlinear  transient  analysis  of  a  rigid  rotor  on  oil  journal 
bearings  under  (1)  a  unidirectional  constant  load,  (2)  a  unidirectional  peri¬ 
odic  load,  and  (3)  variable  rotating  load  Is  performed.  A  number  of  trajecto¬ 
ries  have  been  obtained  with  the  aid  of  a  high-speed  digital  computer  and 
graphi cs . 
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NOMENCLATURE 

radial  clearance 
journal  diameter 
eccentri cl ty 
hydrodynamic  forces 

Ff-  =  FfC^/ncoR^L,  F0  =  F9C2/riuR3L  (dimensionless) 

film  thickness,  h  =  h/C  (dimensionless) 

mass  parameter,  M  =  MCw^/Wq  dimensionless) 

film  pressure,  p  =  pC^/nuR^L 

journal  radius 

dimensionless  time,  T  =  upt 


2 


t 

X  ,  2,  0,  2 

w,  w 

Wo.  Wo 

e 

n 


time 

coordinates,  9  =  x/R,  z  =  z/(L/2)  (dimensionless) 
load,  W  =  WC^/ncjR^L  (dimensionless) 
steady-state  load,  Wq  =  WC^/nuR^L  (dimensionless) 
eccentricity  ratio,  e  =  e/C  (dimensionless) 
absolute  viscosity  of  oil 
attitude  angle 
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Angular  nat'’"  at  which  film  comniences  aim  Laviiates,  lespec- 

tlvely,  from  minimum  film  thickness 
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0]  a  rr  +  f  ]  ,  02  «  ir  +  H)2 
whirl  ratio,  Q  =  oip/oi  (dimensionless) 
angular  velocity  of  journal 
angular  velocity  of  whirl 


THEORY 

The  basic  differential  equation  for  pressure  distribution  In  the  bearing 
clearance  under  dynamic  conditions  can  be  written  as  (See  f1g.  1.) 

8  A  3  3P\  3  /'uS  c  o  f  0  §2^  M  19  Sh  m 

3x  V  3xj  3z  V  dz)  “  (u  -  2  gj  3x  ^  (1 

Equation  ^1)  when  nondimenslor ' jz-?d  with  the  following  substitutions; 

9  =  x/R.  z  =  z/(L/2),  h  =  h/C,  .  -  X^/nuR^,  T  =  oipt,  and  Q  =  uip/oi,  will  read 
as 


le  t)  *  S)  ■  <'  if)  i  ‘  I? 

The  film  thickness  (dimensionless)  Is 

h  =  1  +  e  cos  0 

with  the  use  of  equation  (3),  equation  (2)  can  be  written  as 


(2) 


(3),oT - - 


30 
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=  6(1  -  20cp)(-c  sin  9)  +  12nc  cos  9 


where  <p  =  3<p/3T  and  e  =  3e/3T.  In  this  study  we  have  used  the  Reynolds 
boundary  conditions,  which  are  given  by: 
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vvf.ere  0]  and  02  are  the  angular  coordinates  at  which  the  film  commences  and 
refurms,  respectively. 


Equation  (4)  1s  solved  numerically  for  pressure  by  a  finite  difference 
method  with  a  successive  ove;  r el axation  scheme  satisfying  the  above  boundary 
conditions.  Initially,  e  and  ^  are  set  equal  to  2ero  to  obtain  the  steady- 
state  hydrodynamic  forces  Ff  and  F0. 

The  forces  are  computto  from 
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Releasing  the  journal  from  the  steady-state  position,  one  can  compute  e, 
e,  and  for  the  next  time  step  by  solving  the  following  equations  of 
motion ; 


MC 


^d^e 


2  ■  ^Vdt. 


MC 


Ldt 


=  F^  -I-  W  cos  cp 


=  Fq  -  W  sin  (p 


(7) 


(8) 


The  dimensionless  form  of  equations  (7)  and  (8)  are: 


Mn^e  -  Rfl^e<p^  -  -  cos  !p  =  0 

W. 


(9) 
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(10) 


where 


MQ^e  +  2RQ^e<p  -  ■— -t-  s1n  <p  =  0 


M  inrt  n 

M  =  -j^  and  Wq  = 

o  nwR  L 


The  steady-state  load  Wq  Is  obtained  by  letting  e  and  9  equal  to  zero. 

Equations  (9)  and  (10)  are  second-order  differential  equations  In  e  and 
cp.  These  are  solved  by  using  a  fourth-order  Runge-Kutta  method  for  constant 
values  of  Q,  M,  Ff,  and  F0. 

In  the  following  section  three  types  of  load  are  considered. 


A  Unidirectional  Constant  Load 

Assuming  the  initial  conditions  given  In  table  I,  the  hydrodynamic  forces 
under  steady-state  condition  are  found.  Equations  of  motion  are  solved  to 
obtain  e,  e,  9,  and  9  for  the  subsequent  time  step.  Now  the  new  value  of 
c,  c,  and  9  are  Introduced  In  equations  (3)  and  (4)  to  determine  the  hydro- 
dynamic  forces.  These  forces  along  with  the  steady-state  load,  mass  parameter 
and  whirl  ratio  are  utilized  for  the  solution  of  equations  (9)  and  (10).  The 
process  Is  repeated  until  we  get  a  trajectory  that  describes  the  status  of  the 
system. 


A  Unidirectional  Periodic  Load 
The  applied  load  was  assumed  to  be 

W  =  Wq  [1  +  s1n(l  t)]  (11 ) 

where  Wq  Is  the  applied  steady-state  load  consistent  with  the  eccentricity, 
cq  =  0.8.  At  each  time  ste£  a  new  value  of  W  was  calculated  from  equa¬ 
tion  (11).  For  values  of  M  and  Q  given  In  table  I  and  using  W,  F^,  and 
F9,  equations  (9)  and  (10)  were  solved  for  e,  e,  9,  and  9  for  the  next  time 
step.  The  rest  of  the  procedure  is  repeated  as  described. 


Variable  Rotating  Load 

This  type  of  loading  Is  Interesting  since  It  pertains  to  an  engine  bear¬ 
ing.  The  data  used  for  the  analysis  refer  to  the  Ruston  and  Hornsby  6  VEB-X 
MKIII  connecting  rod  bearing  (ref.  5)  and  Is  listed  In  table  I.  The  polar  load 
diagram  Is  represented  by  figure  2.  For  the  magnitude  of  load  at  different 
positions  of  crank  angle  the  reader  may  refer  to  appendix  4.1  of  the  paper 
(ref.  5).  The  time  step  AT  »  ir/18  was  taken  In  this  case  to  match  with  the 
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given  applied  load  at  10°  crank  angle  Interval.  The  resultant  applied  load  was 
nond  1  men i i vona  1 1  zed  using  the  expression  W  =  WC^/pcjR^L.  This  load  was  Intro 
duced  In  the  equations  of  motion  to  determine  the  position  of  the  journal  cen¬ 
ter  (e,  £,  9,  and  <p>.  The  hydrodynamic  forces  Fp  and  F9  were,  however, 
computed  after  solving  equation  (4). 


RESULTS  AND  DISCUSSION 

To  assure  ourselves  that  our  formulation  Is  consistent  with  reference  2, 
our  results  are  shown  In  figure  3(a)  and  compared  to  those  of  reference  2  shown 
1n  figure  3(b).  The  comparison  Indicates  very  good  agreement  for  the  condi¬ 
tions  stated  in  the  figure.  In  figure  4(a)  a  typical  three-dimensional  pres¬ 
sure  distribution  for  a  particular  time  step  Is  shown.  From  this  figure  one 
can  see  the  cavltated  region  and  the  variation  of  pressure  throughout  the  com¬ 
plete  clearance  space  of  the  bearing.  Figure  5  gives  the  journal  locus  when 
the  bearing  Is  stable  under  the  action  of  constant  unidirectional  load.  In 
figure  3.  the  journal  trajectory  within  the  clearance  circle  Is  shown  for  an 
unstable  bearing.  The  journal  locus  reaches  a  limit  cycle.  This  type  of 
ODservatlon  h?-  been  made  by  many  workers  (refs.  6  to  8)  while  dealing  with 
short  and  infinitely  long  bearings. 

Having  compared  the  result  of  the  present  solution  with  that  of  Akers 
et  al .  (tef.  2),  an  attempt  was  made  to  study  the  effect  of  periodic  and  varia¬ 
ble  load.  Figure  6  shows  the  journal  locus  for  a  periodic  load  superimposed  on 
the  constant  load  given  the  same  as  those  of  a  unidirectional  constant  load. 
From  these  two  figures  it  may  be  seen  that  a  bearing  which  Is  stable  under  the 
action  of  a  unidirectional  constant  load  can  be  made  unstable  when  a  periodic 
load  is  applied.  In  the  latter  case  the  journal  locus  reaches  a  limit  cycle. 

The  journal  center  trajectory  of  the  connecting  rod  bearing  (variable 
rotating  load)  is  shown  in  figure  7.  The  trajectory  Is  very  complex:  It  nei¬ 
ther  tends  to  go  to  an  equilibrium  point  nor  to  reach  a  limit  cycle.  It  may  be 
mentioned  that  the  value  of  the  minimum  film  thickness  for  this  case  exceeds 
that  obtained  by  others  using  the  mobility  method  of  solution.  The  mobility 
method  (ref.  3),  however,  does  not  consider  the  stability  of  the  system. 

The  present  solution  does  not  take  account  of  the  transport  of  fluid 
through  the  cavltated  region  and  consequently  does  not  take  full  account  of  the 
oil  film  history.  A  subsequent  analysis  has  been  undertaken  that  will  Include 
the  oil  film  history  effects  for  these  dynamic  conditions  and  will  be  compared 
at  a  later  date.  A  proper  accounting  of  the  mass  flow  In  time  Is  believed  to 
be  Important  In  dynamically  loaded  bearings  and  should  Influence  the  trajectory 
of  the  rotor  in  a  more  realistic  way. 


CONCLUSIONS 

From  the  nonlinear  transient  analysis  of  an  oil-film  journal  bearing  under 
different  dynamic  loads  with  Reynolds  type  boundary  conditions,  the  following 
conclusions  are  evident. 
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1.  Although  the  dnalysls  1s  costly  In  terms  of  computer  time,  It  gives  the 
orbital  trajectory  within  the  clearance  circle  which  Is  not  obtainable  using  a 
simplified  (linearized)  theory. 

2.  As  it  is  shown  by  others  while  dealing  with  short  and  Infinitely  long 
bearing  theory,  the  journal  locus  for  a  finite  bearing  ends  In  a  limit  cycle 
for  an  unstable  system. 

3.  A  stable  journal  with  a  constant  unidirection  load  can  be  made  unstable 
by  superimposing  a  periodic  load  on  the  system. 

4.  A  more  correct  (cavitated)  boundary  condition  using  the  prehistory  of 
film  may  show  some  interesting  phenomena  which  are  not  revealed  in  the  present 
method  of  solution. 
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TABLE  I.  -  INPUT  DATA  FOR  VARIOUS  CONDITIONS  OF  LOAD 


Cond i t i on 

L/D 

*^0 

M 

Q 

Un i d i recti onal 
constant 

1  .0 

0.8 

5. 

0.5 

Unidi rectional 
period! c 

1  .0 

0.8 

5. 

0.5 

Vari  abl  e 
rotati ng 

0.562 

0.724 

5. 

1.0 

Q  oo(RAD/S)  0(M)  C/R  n(Pa-s) 


1.0  62.84  0.2  0.0008  0.15 


I  !UJRt  It.  A  TYPICAL  PRf  SSURt  DISTRIBUTION  (l/D  1.0,  Eg  =  0.8.  FIGURE  5.  -  JOURNAL  CENTER  TRAJECTORY  FOR  A  UNIDIRECTIONAL  CONSTANT 

10,.  LOAD  <L/D  =  1.0,  E„  =  0.8,  H  =  5.  IT  =  0.5). 


FIGURE  b.  -  JOURNAI  CINTIR  IRAJECTORY  FOR  A  UNIDIRECTIONAL  PERIODIC  FIGURE  /.  -  JOURNAL  CENTER  TRAJECTORY  FOR  VARIABLE  ROTATING  LOAD. 

lOAD  (l/D  1.0,  E_  0.8.  M  S.  II  -  0.5) 
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^  Most  published  work  relating  to  dynamically  loaded  journal  bearings  are  directed  to  determining  the  minimum  film 
thickness  from  the  predicted  journal  trajectories.  These  do  not  give  any  information  about  the  sub;  ynchronous 
whirl  stability  of  journal  bearing  system.*l(^since  they  do  not  consider  the  equations  of  motion.  It  is,  however, 
necessary  to  know  whether  the  bearing  systemVperation  is  stable  or  not  under  such  an  operating  condition.  The 
purpose  of  t  e  present  paper  is  to  analyze  the  stability  characteristics  of  the  system.  A  linearized  j  erturbation 
theory  about  the  equilibrium  point  can  predict  the  threshold  of  stability;  however  it  does  not  indiede  postwhirl 
orbit  detail,  "he  linearized  method  may  indicate  that  a  bearing  is  unstabl^for  a  given  operating  con  lition  whereas 
the  nonlinear  analysis  may  indicate  that  it  forms  a  stable^imit  cycle.  For  this  reason,  a  nonlinear  transient 
analysis  of  a  rigid  rotor  supported  on  oil  journal  bearings  under  (1)  a  unidirectional  constant  load,  (2)  a  unidireC' 
tional  periodic  load,  and  (3)  variable  rotating  load  arc  performed.  In  this  paper,  the  hydr  xlynamic  forces  are 
calculated  after  solving  the  time-dependent  Reynolds  equation  by  a  finite  difference  meth 'd  with  a  successive 
overrelaxation  scheme.  Using  these  forces,  equations  of  motion  arc  solved  by  the  fourth-  rder  Ru;  ge-Kutta 
method  to  prccict  the  transient  behavior  of  the  rotor.  With  the  aid  of  a  high-speed  digital  compute;  and  graphics, 
the  journal  trajectories  arc  obtained  for  several  different  operating  condition, s, 
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